DFF 429/11/05 



The role of orthogonal polynomials in the six-vertex 
model and its combinatorial applications* 



F. Colomo and A. G. Pronko** 

I.N.F.N., Sezione di Firenze, and Dipartimento di Fisica, Universitd di Firenze, 
Via G. Sansone 1, 50019 Sesto Fiorentino (FI), Italy 



Abstract 

The Hankel determinant representations for the partition function and boundary cor- 
relation functions of the six-vertex model with domain wall boundary conditions are in- 
vestigated by the methods of orthogonal polynomial theory. For specific values of the 
parameters of the model, corresponding to 1-, 2- and 3-enumerations of Alternating 
Sign Matrices (ASMs), these polynomials specialize to classical ones (Continuous Hahn, 
Meixner-PoUaczek, and Continuous Dual Hahn, respectively). As a consequence, a uni- 
fied and simplified treatment of ASMs enumerations turns out to be possible, leading also 
to some new results such as the refined 3-enumerations of ASMs. Furthermore, the use 
of orthogonal polynomials allows us to express, for generic values of the parameters of 
the model, the partition function of the (partially) inhomogeneous model in terms of the 
one-point boundary correlation functions of the homogeneous one. 
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1 Introduction 



We shall consider here the six-vertex model on a square lattice with domain wall bound- 
ary conditions (DWBC). The model, in its inhomogeneous formulation, i.e., with position 
dependent Boltzmann weights, was originally proposed in [1], within the theory of corre- 
lation functions of quantum integrable models, in the framework of the quantum inverse 
scattering method [2]. It was subsequently solved in [3], where a determinant formula 
for the partition function was obtained and proven (see also [4]). Analogous determinant 
formulae have been given also for the one- and two-point boundary correlation func- 
tions [5, 6]. In its homogeneous version, the six- vertex model with DWBC admits usual 
interpretation as a model of statistical mechanics with fixed boundary conditions, and it 
may be seen as a variation of the original model with periodic boundary conditions [7, 8]. 
The model is known to be closely related with the problems of enumeration of alternat- 
ing sign matrices (ASMs) and domino tilings (see book [9] for a nice review). It should 
be mentioned that ASM enumerations in turn emerge, via Razumov-Stroganov conjec- 
ture [10, 11], in the context of some quantum spin chains and loop models; for recent 
works, see for instance [12-14] and references therein. 

Till now, specific results for the six-vertex model with DWBC at particular values of 
its parameters (obtained mainly in application to ASM enumerations) were derived from 
general results for the inhomogeneous version, first specializing the parameters to the 
considered case, and next performing the homogeneous limit. Each time, the homogeneous 
limit was an hard task on its own right, and a specific approach was devised to work it out 
in each single case [15-17]. More recently, especially in the context of Razumov-Stroganov 
conjecture, where the model is considered at the so-called ice-point, the homogeneous limit 
of the partition function is often performed only on a subset of the spectral parameters, the 
remaining ones being reinterpreted as the variables of some generating function encoding 
some peculiar properties of the model [17, 18]. 

Here we consider a difi^erent approach, which turns out to be rather convenient: we 
start directly from the homogeneous limit representations, worked out once for all for the 
model with generic vertex weights. Such 'Hankel determinant' representations has been 
derived in [4] for the partition function, and in [5] and [6], respectively, for the one- and 
two-point boundary correlators. Following a rather standard procedure, these quantities 
can then be expressed in terms of orthogonal polynomials. For specific values of the 
vertex weights corresponding to the so-called 'ice-point', 'free-fermion line' and 'dual ice- 
point' of the six-vertex model, these orthogonal polynomials specialize to classical ones 
(i.e., of hypergeometric type, or, equivalently, belonging to the Askey scheme, see for 
example [19]), allowing the evaluation of the partition function and boundary correlators 
in closed form. These three cases correspond to the 1-, 2- and 3-enumerations of ASMs, 
respectively; a unified and simplified treatment of known ASMs is thus provided, leading 
also to some new results such as the refined 3-enumerations of ASMs. 

For generic values of the vertex weights, such orthogonal polynomial representations 
allow to express the two-point boundary correlator in terms of the analogous one-point 
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Figure 1: The six allowed types of vertices in terms of arrows (first row), in terms of lines 
(second row), and their Boltzmann weights (third row). 

boundary correlator. Such relationship can in turn be understood as follows: the partition 
function of the model with two inhomogeneities may be written in terms of the one point 
boundary correlator solely. This simple result can be extended further: indeed here 
we show that it is possible to 'get away' from the homogeneous limit and recover the 
structure of the inhomogeneous partition function, thus generalizing previous formulae of 
Stroganov [17] to a generic number of unfixed spectral parameters. Specifically, we provide 
a simple expression for the partition function of the (partially) inhomogeneous model in 
terms of the one-point boundary correlators of the homogeneous model; our result is valid 
for any choice of the crossing parameter. This somehow simplifies and render explicit 
previous results discussed in [18], where the partition function of the inhomogeneous 
model (at ice-point) was expressed in terms of Shur functions. 

2 The six-vertex model with DWBC 

The six-vertex model, which was originally proposed as a model of two-dimensional ice 
(hence the alternative denomination: 'square ice'), is formulated on a square lattice with 
arrows lying on edges, and obeying the so-called 'ice-rule', namely, the only admitted con- 
figurations are such that there are always two arrows pointing away from, and two arrows 
pointing into, each lattice vertex. An equivalent and graphically simpler description of 
the configurations of the model can be given in terms of lines fiowing through the vertices: 
for each arrow pointing downward or to the left, draw a thick line on the corresponding 
link. The six possible vertex states and the Boltzmann weights Wi assigned to each vertex 
according to its state i {i = 1, . . . , 6) are shown in figure [T] We shall presently restrict 
ourselves to the homogeneous version of the model, where the Boltzmann weights are site 
independent. We shall however turn to the inhomogeneous version of the model in the 
last section. 

The DWBC are imposed on the N x N square lattice by fixing the direction of all 
arrows on the boundaries in a specific way. Namely, the vertical arrows on the top and 
bottom of the lattice point inward, while the horizontal arrows on the left and right sides 
point outward. Equivalently, a generic configuration of the model with DWBC can be 
depicted by lines fiowing from the upper boundary to the left one. This line picture 
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Figure 2: One of the possible configurations of tlie model with DWBC, in the caseiV = 5: 
(a) in terms of arrows; (b) in terms of lines. 



(besides taking into account the 'ice-rule' in an automated way) is intuitively closer to 
ASMs recalled in the next section. A possible state of the model both in terms of arrows 
and of lines is shown in figure |21 

The partition function is defined, as usual, as a sum over all possible arrow con- 
figurations, compatible with the imposed DWBC, each configuration being assigned its 
Boltzmann weight, given as the product of all the corresponding vertex weights. 



arrow configurations i=l 
with DWBC 



Here rii denotes the number of vertices in the state i in each arrow configuration {rti -(- 
■■■ + nf, = N^). 

The six-vertex model with DWBC can be considered, with no loss of generality, with 
its weights invariant under the simultaneous reversal of all arrows, 

Wi = W2 =■ a, Ws = W4 =: b, = Wq =: c. (2.2) 

Under different choices of Boltzmann weights the six-vertex model exhibits different be- 
haviours, according to the value of the parameter A, defined as 

It is well known that there are three physical regions or phases for the six-vertex model: 
the ferroelectric phase, A > 1; the anti-ferroelectric phase, A < —1; and, the disordered 
phase, — 1 < A < 1. Here we restrict ourselves to the disordered phase, where the 
Boltzmann weights are conveniently parameterized as 

a = sin(A + ?7), 6 = sin(A — r^), c = sin2r]. (2-4) 



With this choice one has A = cos 27]. The parameter A is the so-called spectral parameter 
and r) is the crossing parameter. The physical requirement of positive Boltzmann weights. 



in the disordered regime, restricts the values of the crossing and spectral parameters to 
< r/ < 7r/2 and rj < X < n — rj. 

An exact representation for the partition function (for generic weights, even com- 
plex) was obtained in [4]. When the weights are parameterized according to ()2.4p such 
representation reads 

= - V) ^;^(^ + ^ (2.5) 

where ^ is an x iV matrix with entries 

cP,, = ai+V(A,r7), ^iX,v)= . ''^"^1,^ V (2.6) 

sm(A — Tj) sm(A + rj) 

Here and in the following we use the convention that indices of N x N matrices run over 
the values j, k = 0, . . . , N — 1. 

This formula for the partition function has been obtained as the homogeneous limit 
of a more general formula for the inhomogeneous six-vertex model with DWBC The 
inhomogeneous model, with site-dependent weights, is defined by introducing two sets of 
spectral parameters {Xa}a=i and {2^/3)^=1, such that the weights of the vertex lying at 
the intersection of the a-th column with the /3-th row depends on — z/^ rather than 
simply on A, still through formulae ()2.4|) . The inhomogeneous model can be fruitfully 
investigated through the Quantum Inverse Scattering Method, see papers [1, 3, 4, 20] and 
book [2] for details. As a result, the partition function of the inhomogeneous model is 
represented in terms of certain determinant formula which, however, requires some effort in 
the study of its homogeneous limit, — *• and Aq, — ^ A, since in this limit the determinant 
possesses A^^ — A^ zeros that are cancelled by the same number of singularities coming 
from the pre-factor. A recipe for taking such a limit was explained in [4] where formula 
fl2.5|) was originally obtained. Subsequently, formula ()2.5|) was used in papers [21, 22] 
to investigate the thermodynamic limit, A^ — > 00, of the partition function. In these 
studies the Hankel structure of the determinant appearing in fj2.5|) . a natural outcome of 
the homogeneous limit procedure, was exploited through its relation with the Toda chain 
differential equation and with matrix models. 

In addition to the partition function, in the following we shall discuss one- and two- 
point boundary correlation functions as well. In general, two kinds of one-point correlation 
functions can be considered in the six-vertex model: the first one ('polarization') is the 
probability to find an arrow on a given edge in a particular state, while the second one is 
the probability to find a given vertex in some state i. If one restricts to edges or vertices 
adjacent to the boundary, then such correlation functions are called boundary correlation 
functions. Following the notations of paper [5] , where these boundary correlation functions 
were studied, let G^'' denote the probability that an arrow on the first row and between 
the r-th and (r -|- l)-th columns (enumerated from the right) points left (or, in the line 
language, that there is a thick line on this edge), and let H^^ denote the probability that 
the first vertex in the r-th column (counted from the right) is in the state i = 5 (or that 
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the thick hne flows from the top to the left), see Figs. Q and El The first correlation 

(r) 

function, G)^ , is, in fact, the boundary polarization, whose interpretation is more direct 

(r) 

from a physical point of view, while the second one, H)^ , is closely related to the refined 
enumerations of ASMs. It is easy to see that, due to DWBC, the two correlation functions 
are related to each other as follows 

GP = HP + Hi;-'^ + ... + Hii\ (2.7) 

In [5] both correlation functions were computed using Quantum Inverse Scattering 
Method for the inhomogeneous six-vertex model. In the homogeneous limit, which is the 
situation we are interested in here, determinant formulae generalizing ()2.5j) were found for 

(r) 

these correlation functions. For instance, for -ff)^ , the following expression was derived 

^(r) ^ (iV-l)!sin(2r/) det^^ 

[sin(A + r/)]'^[sin(A-r/)]^-'^+' det^^ 

where the matrix differs from the matrix equation ()2.6p . just by the elements of the 
last column 



[sm(£: + A - r])]^-! 



■ (2.9) 

£ = 



A similar expression is valid for G'^^ as well. In what follows we shall focus on h'^^] 

(r) 



the results for G]^ will follow immediately from relation ()2.7j) . From the DWBC it 

N 



immediately follows that = 1, and therefore, from (|2.7p . correlation function H.'f'^ 
satisfies 



Y.eV = \. (2.10) 



r=l 



In the following this normalization condition will be used in application to the generating 



function of defined as 



N 

i/^(n):=5^/7f-'-+V-\ (2.11) 



Let us finally recall the two-point boundary correlation function. Again several differ- 
ent definitions are sensible. We shall focus here, for the sake of simplicity and definiteness, 
on a specific two-point generalization of the one-point correlation function describ- 
ing the probability of finding vertices of type z = 5 on the opposite, top and bottom, 
boundaries. More precisely, we define H^^^ as the probability of finding vertices of type 
2 = 5 both at the ri-th position of the first row and at the r2-th position of the last row 
(ri and r2 are counted from the right). 
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In [6] such correlation function was computed using Quantum Inverse Scattering 
Method for the inhomogeneous six-vertex model. In the homogeneous limit, which is 
the situation we are interested in here, the following determinant formula was derived: 



H 



{ri,r2) 
N 



{N-l)\{N-2)\ sm\2r]) 



[sin(A + ?7)] [sin(A — 7] 

detN [^j,k 



A''+r2— ri+l 



detAT^ 



h 



0<j<N-l,0<k<N-3 



(r-i,r2) 
N 



(^1,^2) 



[2.12) 



£l=£2=0 



Here the function /^^^'^^■'(ei, £2) is defined as follows: 



sin£i)^-'^i[sin(£i - 2r])Y^-\sme2r-''^[sm{e2 + 2r/)] 



lr-2-l 



sin(£2 -ei + 2r/)[sin(ei + A - r7)]^-2[sin(e2 + A + t])]^-^ 
The matrix <P is defined in (j2.6j) . 



(2.13) 



It is obvious from the definitions of H^'^ and H^^''^', that 



r(^i:''2) 
'at ; 



N N 
ri=l r2=l 



N 



H 



(r) 
N ; 



iri,r2) 
N 



1. 



:2.i4) 



ri,r2=l 



In dealing with the two-point boundary function it will be convenient to use the corre- 
sponding generating function 



N 



(Ar-r+l,s) r-l s-l 



U V 



(2.15) 



r=l,s=l 



Note that ()2.14j) simply implies Hn{1,v) = Hn{v, 1) = Hjsiiy) where Hn{v) is defined by 
I^J^; we also have Hn{1, 1) = 1. 



3 Orthogonal polynomial representation 
3.1 Preliminaries 

In this section the previously reviewed determinant representations for the partition func- 
tion and the one- and two-point boundary correlation functions will be analysed by making 
use of the orthogonal polynomial theory, along the lines proposed in paper [6, 23]. We 
start with recalling some very standard and well-known facts from the general theory. 
Let {P„(a;)}^o be a set of polynomials, with non- vanishing leading coefficient 

P„(x) = + . . . , K„ ^ 0, (3.1) 

and orthogonal on the real axis with respect to some weight /i(x), 

POO 

/ Pn{x) Pm{x) Ai(x) dx = kjnm ■ (3.2) 
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Let Cn denote n-th moment of the weight /i(a;), i.e. 

x"/i(x)dx, n = 0, 1, . . . 



(3.3) 



and let us consider the (n + 1) x (n + 1) determinant 





Co 


Cl 




Dn = 


Cl 


C2 . . 






C-n 


Cn+1 • • 


■ C2n 



(3.4) 



Using the orthogonahty condition ()3.2|1 and well-known properties of determinants, one 
can easily derive the following formula 



Dr 



k=0 1^ 



(3.5) 



This formula can be used for computation of Hankel determinants, provided the orthogo- 
nal polynomials {P„(a;)}^Q are known. On the other hand, the polynomials {-Pn(2^)}^o 
can in turn be expressed as determinants. For later use let us introduce the notation 





Co 


Cl 


■ ■ ^n—k 1 1 


. 1 






Cl 


C2 


C„_fc_|_l Xi X2 ■ 


• Xk 


(3.6) 






Cn+1 


■ ■ C2n—k ^2 


■ ^k 




so that Dn'' = Dn, and D^}-^^{xi 


• • • 




A{xi, . . . , Xn), where A{xi, . 


. . ,Xn) denotes 



the Vandermonde determinant of n variables. For the polynomials one can find that 



Pn{x) 



D 



(3.7) 



n-l 



For a proof, see for example, [29]. 

Relation ()3.7p can be read off inversely thus giving an expression for the determinant 
Dn\x) in terms of the polynomials Pn{x). Taking into account that (see ()3.5|) ) 

h D 

"'n -"-^n 

'^^D~ 



n-l 



we can write 



D^r!Hx) 



^ Pnix). 



^ ^ (3.9) 

^n "'n 

Consider now the case of Di^) (a;i, X2). It is clear that the term of the highest powers 
on both xi and X2 is just Dn-2{x2X^~^ — 0:^X2"^); following [29], it can be shown that 

Dn-2 



D^i 



Xl,X2 



^n^^n—l 



[Pn^l{Xi)Pn{x2) - Pn{Xl)Pn-l{x2)] 



(3.10) 



Again using (jH.Hj) . we write 



Dr 



hnhn—1 
"'n"'7i— 1 



[Pn-l{Xi)Pn{x2) - P„ (Xi )P„_i (Xs)] 



P„_l(Xi) P„-l(x2) 
PniXl) Pn{x2) 



(3.11: 



This formula can be easily extended to the general case of Dn\xi, . . . ,Xk)- We shall 
now use all these formulae to provide completely general 'orthogonal polynomial repre- 
sentations' for the partition function and the one- and two-point boundary correlation 
functions. 



3.2 The partition function 

First we note, following papers [22, 23], that the matrix <P entering the expressions for 
the homogenous model partition function and the boundary correlation functions can be 
related with orthogonal polynomials using the integral representation 

''^[^'^\^ r = r e-(^-/^) dx. (3.12) 

sin(A — r]) sin(A + r]) sinh(7rx/2) 

This formula is valid ifO<77<7r/2 and 1] < X < n — r]] these values of A and 1] correspond 
to the so-called disordered regime of the six-vertex model. When considering other regime 
of the model, the measure of the corresponding polynomials becomes discrete, see [22], 
but the present procedure may nevertheless be considered, modulo trivial modifications. 
Indeed, it can be easily seen that our results below do not depend on the particular choice 
of the regime, and can be extended to other regimes simply using the proper analytical 
continuation in the parameters A and rj. 

Formula ()3.12j) implies that we have to deal with the set of polynomials which are 
orthogonal with respect to the following weight function 

^^{x)=f,{x■,X,v) = e-^'--/'^^^^^. (3.13) 

smh(7ra;/2) 

The corresponding polynomials Pn(x) = Pn(x; \,ri) also depend on A and rj which are to 
be considered as parameters. In what follows we shall often omit the dependence on A 
and r] where possible. Let us mention the following useful property of these polynomials 

P„ix; A, r]) = (-1)" Pni-x; vr - A, 77). (3.14) 

This property can be easily established in virtue of formula (j3.7|) . It is to be mentioned 
also that both the leading coefficient /€„ = Kn{X,i]) and the normalization constant /i„ = 
hn{X, Tj) are invariant under the substitution A ^ tt — A. The transformation A ^ tt — A is 
related to the crossing symmetry of the six-vertex model, and the previous property will 



9 



have useful consequences in the discussion of the one- and two-point boundary correlation 
functions, discussed below. 

Let us focus now on the partition function. The differential operator in the Hankel 
matrix entering representation ()2.5|) pulls down powers of the integration variable, and 
the matrix itself is immediately rewritten as the matrix of moments of the measure ()3.13|) . 
Due to ()3.5|1 . we readily get 

/ . r, N-l , 

\ ^ J „=o y^-r^n 

where (p = (p{X,r)) is exactly the function defining entries of the matrix see ()2.6p . 

We just want to conclude with a simple remark: whenever the values of the parameters 
A and rj are such that the corresponding orthogonal polynomials Pn{x]X,r]) happen to 
belong to the Askey scheme, the evaluation of the partition function in closed form reduces 
to an elementary computation. As it will be shown in the following, this simple observation 
allows for a direct and straightforward evaluation of ASM enumerations. 



3.3 The one-point boundary correlation function 

Let us now turn to the one-point boundary correlator. Starting from its determinant 
representation, formula (j2.8p . and recalling relation (|3.9|) . we readily write: 



(A^-1)! sin(2r/) Kr,.i{X,v) 



[sin(A + r7)]'[sin(A - r])]^^"^^ hN_,{X,ri) 

(sine)^-''[sin(£-2r/)]' 



X PN_i{ds;X,r])- 



(3.16) 



£ = 



[sin(e + A - r])]^-^ 

This representation does not seem particularly appealing at this level, but it is worth 
noticing that, when the values of A, rj are such that P/v--i(a;, A, 77) reduces to a classical 
polynomial, the previous correlation function may be evaluated exactly, in closed form. 
The use of the properties of the polynomial entering the representation is a crucial in- 
gredient in performing such computation. This will be shown in detail later, but now we 
are interested instead in obtaining another equivalent representation for H^p , by making 
use of the crossing symmetry of the six-vertex model. The use of these two equivalent 
representations for the one-point boundary correlation function will then allow us to 
express the two-point boundary correlation function if^^'^^^ in terms of one-point bound- 
ary correlation functions. 

We recall that the crossing symmetry is the symmetry of the vertex weights under 
reflection with respect to the vertical axis, and simultaneous interchange of the functions 
a and b, which is equivalent to setting A vr — A. Since the lattice with DWBC is invariant 
under the reflection with respect to the vertical axis (modulo reversal of all arrows on the 
horizontal edges), the crossing symmetry implies that the following relation holds 

Hi;\X,v) = Hi,''-^^'\n-X,r^). (3.17) 



10 



Consider expression ()2.8|) for the one-point function H^^\ Due to ()3.9|) we can rewrite 
it as 

^(,) _ (AT-l)! sin(2r7) fi:jv-i(A,r/) 



[ sin(A + 7])Y[ sin( A - r/)] /^Jv-i (A, r/) 

(sine)^-^'[sin(£-2r7)]^- 



X PN-l{de] X,r])- 



[sin(e + A - r/)] 



N-l 



. (3.18) 



£ = 



Taking into account ()3.14p and the properties of the leading coefficient Kn{X,ri) and the 
normahzation constant hn{X,ri) mentioned above, it can be easily seen that from ()3.16|) 
and (|3.17|) the following expression is valid as well 



rWn-^N- (A^ - 1)! sin(2?7) KN-iiX,r]) 



[ sin(A + v)Y[ sm( A - r])] /^iv-i (A, v) 

{smeY-^[sm{e + 2ri)y 



[sin{e + X + r])] 



N-l 



(3.19) 



£ = 



Note that this expression means simply that the limit £ — > in ()3.1(')j) can be changed 
into e ^ 21] without altering the result. 

Now, these two equivalent representations, ()3.16|) and ()3.19|) . can be used in the study 
of the two-point correlation function H^^'^^^ given by expression ()2.12j) . which certainly 
involves similar structures. Before turning to this analysis, let us put the above formulae 
for the one-point function in a more compact and convenient notations. 

We define the functions 

sin(A + r/) sine sin(A - r^) sin(e - 2r/) 

^{^) = — ^^ 7 — 7 Q{^) = 7 — ^ 7 ^ 7) (3. 20) 

sm(A — T]) sm[e — 2rj) sm(2?7) sm(£ + X — r}) 

which are related to each other as 



Q{e) = (3.21) 
uj[e) - 1 



Also we define 



^ sin(A-r7) sine _ sin(A + Vj) sin(g + 2ry) _ 

sin(A + 7]) sin(e + 27]) ' sin{2'q) sin(e + A -|- 77) ' ^ ' ' 

which are in turn related to each other as 

Note, that the functions with tildes are introduced such that 

Lo{e]X,vi) =uj{-e;-n - X.ri), g{e; X,7]) = -g{-e;Tc - X,r]) (3.24) 
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in accordance with the crossing symmetry considerations made above. Additionally, let 
us denote 

Kj,_^{x) = (AT - l)!^^^,P^_i(x) (3.25) 

where (f = (f{X,ri) is exactly the function defining entries of the matrix see fl2.fj|l . In 
these notations formulae ()3.16|) and ()3.19|) for the correlation function H^'' read 



and 



= K^^M)He)] [Q{e)f'' (3.26) 

e=0 



respectively. 



<^ = K^.^{^,)[Co{e)Y-'[g{e)f-' , (3.27) 

£=0 



3.4 The two-point boundary correlation function 

Let us now consider the two-point correlation function H^^'^^\ which is given by formula 
()2.12p . Obviously, function ()2.13|) contains all the structures introduced above apart from 
the factor sin(£2 — + 2?]) standing in the denominator there. However, using the identity 

sin(2?7) sin(£:2 ~ £i + 2//) = sinei sin £2 — sin(£:i — 2//) sin(£:2 + 27]) (3.28) 
it can be easily seen that 

sin(£:i + A — r]) sin(£:2 + A + r/) 1 1 



sm{£2 - El + 2r]) ^0i£i)Q{e2) uji£i)uj{e2) - l 



(3.29) 



Thus, taking into account formula (j3.11|) we can write the two-point correlation function 
in the form 



£2) 



^ [c.(.0]^--[g(.,)]^-^[^(.2)]^-"[g>2)]^-^ ^ ^3_3Q^ 



uj[ei)uj[e2) 



Ei=0,e2=0 



Taking into account that uj{e),Cj{e) — > as e — > we can expand the denominator in 
(I3.3(J|) in power series and it can be easily seen that only the first few terms (actually not 
more than A^) of this expansion will contribute. As a result, in virtue of relations (j3.2H) 
and ()3.23|) . we arrive to the following expression in terms of the one-point functions 

N 

TT{ri,r2) _ / TT(ri-j+l) TTiN-n+j) TT{ri~j) tt{N -r2+j) 

- + Hi^:^E'{;-^-^'^\ (3.31) 
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(r) 

where it is assumed that if r < or r > + 1 then = by definition. 

The last expression may be equivalently expressed in terms of the generating functions, 
defined in (HTI^ . (1^^) : 



Hn(u,v) 



1)Hn{u) ■ vHn-i{v) - uHm-i{u) ■ {v - 1)Hn{v) 



V — u 



uHn-i{u) 
[u - 1)Hn(u) 



u — V 

vHn-i{v) 



(3.32) 



This formula generalize to arbitrary values of the vertex weights the result of paper [17] 
where an equivalent expression was derived in the case A = tt /2 and t] = tt/G, i.e., when 
a = b = c (the so-called ice-point). 



4 A combinatorial application 

4.1 Alternating Sign Matrices enumerations 

An alternating sign matrix (ASM) is a matrix of I's, O's and — I's such that in each row 
and in each column (i) all nonzero entries alternate in sign, and (ii) the first and the last 
nonzero entries are 1. An example of such matrix is 



(4.1) 



There are many nice results concerning ASMs, for a review, see book [9]. Many of 
these results have been first formulated as conjectures which were subsequently proved 
by different methods. 

The most celebrated result concerns the total number A{N) oi N x N ASMs. It was 
conjectured in papers [24, 25] and proved in papers [15, 26] that 



/o 


1 








o\ 


1 


-1 


1 

















1 








1 











^0 











V 



N 



A{N) = n 



(3A;-2)!(A;- 1) 



N 

n 



(3A; - 2) 



L(2A:-l)!(2fc-2)! l\i2N-k) 



(4.2) 



A possible generalization of the previous problem consists in considering weighted 
enumerations, or x-enumerations of ASMs. In the x-enumeration matrices are counted 
with a weight where k is the total number of ' — 1' entries in a matrix. The number 
of x-enumerated ASMs is denoted traditionally as A{N;x). The extension of the x = 1 
result above to the case of generic x is not known, but for a few nontrivial cases, namely 
X = 2 and x = 3 [15, 24, 25], closed expressions for x-enumerations are known (note that 
the case x = is trivial, since assigning a vanishing weight to each ' — 1' entry restricts 
enumeration to the sole permutation matrices: A{n] 0) = n\). 
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+ + + + + + 

1 -1 

Figure 3: Vertex states — ASM's entries correspondence. 



A further generalization of the previous problems consists in the so-called refined 
enumerations of ASMs, where one counts the number of N x N ASMs with their sole '1' 
of the last column at the r-th entry. The refined enumeration can be naturally extended 
to be also an x-enumeration. The standard notation for the refined x-enumeration is 
A{N,r;x); in the case x = 1 one writes simply A{N,r) just like A{N) for the total 
number of ASMs. The answer for the refined x-enumeration is known only for the two 
values X = 1 [16, 24, 25] and x = 2 [25, 27]. 

Such enumerations can be further generalized to the doubly refined weighted countings, 
A{N, r, s] x) where one counts the number of N x N ASMs with their sole '1' of the first 
and of the last column at the r-th and {N — s + l)-th entries, respectively. An answer for 
A[N, r, s] x) in the case x = 1 was found in [17]. 

The most direct way to derive ASM enumerations is based on the one-to-one corre- 
spondence between N x N ASMs and configurations of the six-vertex model on N x N 
lattice with DWBC, which has been pointed out in [27, 28], and applied for the first time 
in [15]. The correspondence between matrix entries and vertices is depicted in figure 01 
For example, matrix 1)4.11) corresponds to the configuration of figure |21 and vice versa. 

As an immediate consequence of this correspondence, ASM enumeration is exactly 
given by the partition function of square ice, when all vertex weights are set equal to 
unity. More generally, the number of '—1' entries in a given ASM being equal to the 
number of vertices of type 6 (see figure Ej), and the number of vertex of type 5 and 6 being 
constrained by the condition — uq = N, we readily get 

A{N;x) = (1 -x/4)-^V2^-A^/2^^|^^^^^ _ (4^3) 

»7=arcsin(-y^/2) 

Therefore, x-enumeration of ASM corresponds to the computation of the partition func- 
tion of square ice on the subset of parameters space given by a = 6. In this correspondence, 
values of x belonging to the interval (0, 4) corresponds to the disordered regime of the 
model, -1 < A < 1. 

This nice correspondence can be further extended to the refined x-enumeration of 
ASMs. In the language of square ice, the ratio A{N,r;x)/A{N;x) can be rephrased as 
the probability of finding the unique vertex of type 5 on the first row at the (A^ — r + 1)- 
th site, which is exactly the definition of the boundary correlation function Hj^ 
Explicitly, one has 



A{N,r;x) rr{N-r+l) 



A{N;x) 



A=./2 • (4.4) 

r7=arcsin(Y^/2) 
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Note that for the particular value A = 7r/2, which is the one of interest in ASMs enumer- 
ations, the one-point boundary correlator enjoys the symmetry i7(;) = i7f-^+^). 

Analogously, the doubly refined x-enumeration of ASMs is simply related to the two- 

N 



point boundary correlation function i?^'**^ as follows: 



A{N, r, S; x) _ rN~r+l,s) 

A{N;x) 



A=./2 • (4.5) 

r]=arcsin(y/x / 2) 



It is evident from the previous formulae that 1-, 2- and 3-enumerations of ASMs cor- 
respond to the values rj = n/G, 7r/4, 7r/3, respectively, in the six- vertex model. In these 
three cases, the orthogonal polynomials appearing in the representations of section El for 
the partition function and boundary correlators turn out to specialize to so-called classi- 
cal ones, i.e. appearing into the Askey scheme of hypergeometric orthogonal polynomials. 
Such polynomials can be expressed as terminating hypergeometric series, and hence, their 
characteristic properties, such as orthogonality condition, three-term relation, etc., can 
be worked out explicitely (for full details, and g-analog extensions, see [19]). This hyper- 
geometric structure allows one to evaluate ASM enumerations in a simple and explicit 
way [23]. 



4.2 1-enumerations of ASMs 

To illustrate the method we shall here restrict ourselves to the case of 1-enumerations of 
ASMs. This corresponds to the so-called ice-point, or A = 1/2 symmetric point, of the 
six-vertex model. The parameters of the model are specialized to rj = it /6 and A = 7r/2, 
and correspondingly, the Boltzmann weights assume the values a = b = c = 1. The 
orthogonality weight reads 



/i(x) 



sinh fx 

o 

sinh |x 



4vr2 



1 X 
- + 1- 
3 6 



2 X 
3+^6 



(4.6) 



Direct inspection of some tables of classical orthogonal polynomials, such as [19], allows 
to recognize here a particular specialization of the orthogonality weight for Continuous 
Hahn polynomials, which are defined as 



, , .„(a + c)n(a + d)n ^f-n,n + a + b + c + d-l,a + ix 

Pn{x;a,b,c,d) =1 j 3i^2 , , , 

n! \ a + c , a + a 

and satisfy the orthogonality relation 



(4.7) 



1 

2^ 



Pn{x] a, b, c, d)pm{x] a, b, c, d) T{a + ix) T(b + ix) T{c — ix) T{d — ix) dx 

) 

T{n + a + c) T{n + a + d) V{n + b + c) V{n + b + d) 



{2n + a + b + c + d-l)T{n + a + b + c + d-l)n\ 



Srirr 



(4. 
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Orthogonality condition ()4.8|1 is valid if the parameters a, b, c, d satisfy Re (a, b, c, d) > 0, 
a = c and b = d. Comparing ()4.6|) with ()4.8|) we are naturally led to the choice of 
parameters a = c = 1/3 and b = d = 2/3. Hence, the appropriate polynomials to be 
associated to the Hankel determinant in this case are 



PJx) 



Pn 



x 1 2 1 2 
6' 3'3'3'3 



i"(2/3)„ sF2 



—n, n + 1, l/3 + ix/6 
2/3, 1 



(4.9) 



The normalization constant and the leading coefficient are readily computed: 

2(3n+l)! (2n)! 



hr, 



(2n + 1) 33'^+V2 n! 



6" {n\f 



(4.10) 



Substituting the obtained values of hn and Kn in expression (j3.15|) for the partition 
function, and cancelling whatever possible, we arrive to the following value for the ice 
point partition function 



X=n/2 
ri=n/6 



2 



Ar2 AT-i 

n 

n=0 



(3n + l)!n! 
(2n)! (2n + 1)!' 



(4.11) 



The product expression here gives exactly the total number of ASMs, A{N), since by 
formula (j4.3|l the first factor relates A{N) with the partition function. 



A{N) = (3/4)-^'/2 



A=7r/2 • 
ri=Tr/6 



(4.12) 



Thus, we have easily recovered the celebrated result, Eqn. ()4.2j) . for 1-enumeration of 
ASMs. It is worth noting that the proof presented here is considerably simpler in com- 
parison to those of papers [15, 26]. 

The refined enumeration A{N,r) can be obtained within the same framework, the 
derivation being only slightly more involved. The key ingredient is now that the orthogonal 
polynomials belonging to the Askey scheme are known to satisfy some differential or finite 
difference equation in their variable. In the present case, such equation reads 











I- 


"6 j 


1- 


"6 j 



Pjv-i(a; + 6i) + 



4 

18 9 ^ ^ 



+ 



1 ix 
3^6" 



2 ix 
3^6" 



PN-iix - 6i) = . (4.13) 



First we recall that the refined enumeration is, modulo an obvious overall normaliza- 
tion factor, see ()4.4|1 . nothing but the one-point boundary correlator, which in the present 
case (?7 = 7r/6 and A = 7r/2) reads: 



H 



(r) 
N 



H 



(N-r+l) 
N 



const X \P^^,{deMe)f~'[uj{e)] 



ir-l 



e=0 



(4.14) 
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where Lj{e) and g{6) are given by 

sine sin(e-7r/3) /.1r^ 

^ ^ = —( 7^' ^? ^ = , /o ■ 4.15 

sin(£: — 7r/3) sin(£: + 7r/3) 

Then, exploiting the simple relation P^_i{de ± 6i) = e^^'^P/v-i(<9e)e^^'^, it is easy to see 
that finite difference equation ()4.13|1 implies the condition 

Pjv-i(<9e) [sinSeS^ sin3e + sin2 3e-9A^(A^- 1)] [Q{t)f-^[uj{t)X'^T{t) =0 (4.16) 

where T{e) is at this stage an arbitrary function. Our aim now is to determine the form 
of r(e) in such a way that the last equation can be turned into a recurrence relation in 
r for the boundary correlator h'P . A constructive procedure has been devised in [23] to 
perform this task in general. The derivation is given there in full detail. In the present 
case the resulting recurrence relation reads 

r(r -2N + 1)H^n'^^^ -{r~N){N + r- l)H^j^^ = 0. (4.17) 

This recurrence can be easily solved modulo a normalization constant 

(r) (N + r ~2)\(2N -l-r)\ 

hV = const X (^_//(^_^), ^ (4.18) 

A possible way to fit the normalization condition ()2.1()j) . is to consider the generating 
function H^{z) defined via Eqn. ()2.11|) . The result reads 



(2iV-l)!(2iV-2)! /l-iV,Ar 
(Ar-l)!(3Ar-2)! 2-2N 



(4.19) 



where the proper normalization is easily determined through Chu-Vandermonde identity 

-m , b 

c 

Inspecting the coefficient of z^~^ in ()4.19|) one obtains 

(N+r-2\ (2N-l-r\ 
r(r) _ I Af -1 ) [ N ~l ) 

{ N-1 ) 



^j_(c6)m. (a)^ :=a(a + l)---(a + m- 1). (4.20) 



The refined 1-enumeration of ASMs, conjectured in [25], and first proven in [16], imme- 
diately follows from relation ()4.4j) . 

Let us conclude with the doubly refined 1-enumeration of ASMs, which is readily 
obtained by substituting expression 1)4.211) for the refined enumeration into 1)3.31)) . Anal- 
ogously, the corresponding generating function, Hn{u,v), is immediately obtained by 
substituting ()4.19j) into ()3.32|) . hence recovering the result of paper [17] where a different 
approach was considered. 
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4.3 Other ASMs enumerations 



We have just seen that the present approach allows a very simple, straightforward and 
unified derivation of all answers concerning 1-enumerations of ASMs. These answers were 
already known, but the method exposed above can be extended to all x-enumerations 
where the underlying orthogonal polynomials appear to belong to the Askey scheme. 
Indeed, when this is the case, /i^, i^n are explicitly known, and moreover the considered 
orthogonal polynomials happen to satisfy a finite difference equation in their variable; has 
illustrated above, this is sufficient to compute explicitly the ordinary, refined, and doubly 
refined enumerations. As a matter of fact, it appears that the underlying orthogonal 
polynomials belong to the Askey scheme only for x = 1, 2, 3 (for completeness we mention 
also the trivial case a; = 0, where Meixner-Pollaczek orthogonality measure naturally 
emerges when the = limit is performed in ()3.13|) . and the highly non-trivial case x = 4, 
where the corresponding orthogonal polynomials turns out to be of Bannai-Ito type [30], 
and therefore do not fall within Askey scheme). Our approach allows to recover in a 
straightforward way the results concerning 2-enumerations, and to derive the previously 
unknown answers for the refined and doubly refined 3-enumerations. 

The 2-enumerations correspond to the so-called free-fermion point of the square ice, 
where the parameters of the model assume the values A = 7r/2, r7 = 7r/4. As a matter 
of fact, the whole free-fermion line of the model (A = 0) can be treated at once, with 
the crossing parameter set to the value t] = tt/A, while A is free to vary in the interval 
7r/4 < A < 37r/4. In this case the associated orthogonal polynomials can be recognized as 
a particular specialization of Meixner-Pollaczek polynomials [19]. The computation of the 
partition function and boundary correlators, and thus of the various ASM 2-enumerations, 
can be performed along the lines of the proposed approach [23], and is particularly simple. 
The various answers in this case have been known for some time [4, 5, 25, 27]. 

Let us now discuss the 3-enumerations, which correspond to the so-called dual ice-point 
(A = —1/2) of the six-vertex model, with A = 7r/2 and 7] = ir/S. In this case one can 
recognize that the measure is associated to the Continuous Dual Hahn polynomials [19]. 
A complication however occurs: on one hand these are polynomials in x^ rather than in x, 
and on the other hand the support of the integration measure is restricted to the positive 
half-axis. These inconveniences are however easily circumvented, provided that the cases 
of even and odd N {N being the size of the lattice, or of the ASMs) are treated separately. 
Indeed the set of polynomials Pn(x) associated to the determinant representation for the 
partition function or the boundary correlation functions should be specified differently 
for even or odd n, each of these two cases corresponding to a slightly different choice of 
parameters in the Continuous Dual Hahn polynomials. This is in fact a merely technical 
complication and the general procedure previously outlined may again be applied without 
modification. This has been done in full detail in [23]. We shall here just recall the 
main results. The known answer for the partition function, or, equivalently, for the 
ordinary 3-enumeration of ASMs [15, 25], is readily recovered. Moreover the previously 
unknown answer for the one-point boundary correlator (or, equivalently, for the refined 
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3-enumeration of ASMs) can be worked out, and reads: 



(r) _ B{m,r - 1) + B{m,r - 2) 

(^) _ 2S(m,r-l) + 55(m,r-2) + 2fi(m,r-3) ^' ' 

where the quantities B{m, n), obeying B{m, n) = B{m, 2m — n), are given by 

[n/2] 



B{m, n) 



(2m + l)!m! 
3™ (3m + 2)! 

^=max{0, n—m) 



2m + £- n+l\/'m + £ + l 
m+1 /V m+1 



2"-2^ (4.23) 



for n = 0, 1, . . . , 2m, while they are assumed to vanish otherwise. It should be mentioned 
that this result has first been obtained in [31] through the solution of a second order 
differential equation associated to the generating function Hj^{u), see (12.111) . Such dif- 
ferential equation, which on one hand is equivalent to the recurrence relation which can 
be derived from the orthogonal polynomial approach (i.e. the analogue of (j4.17p for the 
present case), is on the other hand closely related to the functional equation proposed 
in [32], within a different approach to the problem. 

The two-point boundary correlator is immediately obtained simply by inserting the 
previous expression into (|3.31|) . 



5 An alternative representation for the (partially) in- 
homogeneous partition function 

We have till now restricted ourselves to the homogeneous version of the six-vertex model. 
This has allowed us to apply very standard tools from the theory of orthogonal polynomials 
to derive in a simplified and unified way several nontrivial results. Of course, in doing 
so, we have lost the dependence in, let us say, the partition function, on the spectral 
parameters of the model. This dependence implicitly encodes information about the 
correlation functions of the model, and is thus of great interest. We want here to show 
how the results of sectionElfor the boundary correlators, based on the proposed orthogonal 
polynomial representation, can be used to recover the dependence of the partition function 
on the spectral parameters of the (partially) inhomogeneous version of the six-vertex 
model. 

As briefly explained in section |2l when considering the inhomogeneous version of the 
six- vertex model, one associates a spectral parameter (a = 1, . . . , A^) to each row 
of the lattice, and z/^ (/3 = 1, . . . , A^) to each column. We want here to consider a 
particular situation in which all horizontal inhomogeneities 1^13 are set to zero, and a 
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subset of the vertical one is set to some value A. For simplicity, and comparison with 
previous works [17, 18] we shall fix this value to A = 7r/2 (note however that our results 
can be readily extended to the case of generic values of A) . On the other hand 1] will be 
considered generic. We reparametrize the Boltzmann weights ()2.4|) as follows: 



a{z) = zq + 1 , b{z) = z + q . 



c{z) 



z 



where 



^2j(Aa-7r/2) 



(5.1) 



(5.2) 



We have here ignored the parameters which have all been set to zero; the definition 
of the Boltzmann weights in ()5.1|) moreover differs from ()2.4j) by a common overall fac- 
tor. This choice turns out to be more convenient for what follows, and can be absorbed 
in the normalization of the partition function, which is now a function of N variables: 
Zn{zi, . . .,zn). 

In the case when only one inhomogeneity is present, in the first column, the partition 
function can be expressed in terms of the boundary correlator H^^ as follows: 



[z, 1, 



N 

c{z) f a{z) 



b{z) 
6(1) 



r-l 



a[z] 



c(l) \a{l 



H 



(r) 



c(l) \a{l) 



N-l 



HN{h{z)/a{z)) 



(5.3) 



where the generating function Hn{w), see fl2.11|l . is a polynomial of order N — 1 in w. 
Let us now introduce the variable m, related to z as follows 



u 



h{z) 



z + q 



u — q 



a{z) 1 + zq 1 — uq 

and consider the partition function as a function of the variables ui. 

N 



(5.4) 

. , Mn, by introducing 



Zn(ui, 



,Un) 



Zn{zi, 



zn) 



^7V(1, 



n 

a=l 



c(l) fail) 



c{za) \a{za) 



N~l 



(5.5) 



Note that the 'partition function' Zn{ui, 
1. Relation ()5.3|) then reads 

Zn {u, 1, • 



,un) is normalized such that Zn{1, . . . , 1) 



Hn{u) 



(5.6) 



In the same way, following [17, 18], we may rewrite the partition function with two 
inhomogeneities as the generating function ()2.15|) of the two-point boundary correlator. 
Since the latter is in turn expressible in terms one-point boundary correlators, see ()3.32|) . 
we immediately have 



Zn iui,U2, 1, 



1) 



A{ui,U2) 



UiHm-i{Ui) 

[ui - 1)Hn{ui] 



{U2 - 1)Hn{u2) 



(5.7) 
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It is now quite natural to guess that 



ul~'^{Uk - l)HN-k+2{Uk) 

(5.8) 

We note that the right hand side is by construction a symmetric polynomial in mi, . . . , m^, 
of order — 1 in each variable, as it should. Moreover it is evident that the homogeneous 
limit can be performed one variable at a time, say Uk —>■ 1, maintaining the proposed 
structure: formula ()5.8p obviously satisfies 

lim Zn . . . , Uk-i, Uk,l, ■ ■ ■ A) = Zn {ui, . . . , Uk-i, 1,1,...,!) . (5.9) 

It is also straightforward to verify that 

lim Zn (mi, . . . , Wfc-i, Wfc, 1, . . . , 1) = Hn{0) Z^-i (mi, . . . , Wfc-i, !,...,!) , (5.10) 

which is nothing but Korepin recursion relation [1] specialized to the present situation. 
It should however be mentioned that in the present case, where all 'horizontal' spectral 
parameters vp have been set to same value, the uniqueness of the solution of such a 
recursion relation is not guaranteed. A complete derivation of representation ()5.8p can 
however be given by exploiting a set of identities relating the derivatives of function y3(A, 77) 
to the polynomials Hn{u); this will be done elsewhere. 

We would like to conclude by emphasizing that representation (j5.8p constitutes a rather 
wide generalization of previous expressions given in [17, 18]. First of all it generalizes the 
previous formulae, holding in the case of two spectral parameters, to a larger number of 
variables. Moreover, the proposed representation holds for any values of 77, (and can be 
extended straightforwardly to any value of A 7^ 7r/2, by a suitable redefinition of vari- 
able u), expressing in general the partition function in terms of just one-point boundary 
correlators of the corresponding homogeneous model. 
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